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Summary: We present new time-resolved Ha, Call H k K and Mgll h k k spectra of the 
rapidly-rotating KO dwarf star AB Doradus (= HD 36705). The transient absorption features 
seen in the Ha line are also present in the Ca II and Mg II resonance lines. New techniques 
are developed for measuring the average strength of the line absorption along lines of sight 
intersecting the cloud. These techniques also give a measure of the projected cloud area. 
The strength of the resonance-line absorption provides useful new constraints on the column 
densities, projected surface areas, temperatures and internal turbulent velocity dispersions of 

the circumstellar clouds producing the absorption features. At any given time the star appears 

T-o r L 

to be surrounded by at least 6 to 10 clouds with masses in the range 2 4- 6xl0 17 g. The 

T P _ 

clouds appear to have turbulent internal velocity dispersions of order 3*20 comparable 

with the random velocities of discrete filamentary structures in solar quiescent prominences. 
Night-to-night changes in the amount of Ca II resonance line absorption can be explained by 
changes in the amplitude of turbulent motions in the clouds. The corresponding changes in the 

r -- v a ’«'■ f J 

total energy of the internal motions are of order l(r 9 erg per cloud. Changes of this magnitude 
could easily be activated by the frequent energetic (~10 34 erg) X-ray flares seen on this star. 
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1 Introduction 


AB Doradus is a young, spotted KO dwarf with an unusually short axial rotation period of 0.514 
day (Pakull 1981) . It is a flaring coronal X-ray source (Collier Cameron et al 1988; Vilhu & 
Linsky 1987) . Optical and UV spectra show strong chromospheric and transition region line 
emission (Vilhu et al. 1987; Rucinski 1985), indicating high levels of magnetic activity in the 
star s outer atmosphere . Innis et al. (1986) suggested on the basis of its space motion that it 
is a likely member of the Pleiades Moving Group (also known as the Local Association: Eggen 
1973, 1975, 1983a, b). This implies that AB Dor is similar in age to the Pleiades cluster. This 
is supported by its high lithium abundance (Rucinski 1982) and the existence of a population 
of similarly rapidly-rotating and magnetically-active K dwarfs in the Pleiades cluster (van 
Leeuwen et al 1983, 1987; Stauffer et al 1984, 1987). 

In earlier papers in this series we discussed the properties of a system of prominence- 
like condensations of mainly neutral material, located between 3 and 10 stellar radii from the 
rotation axis of AB Dor and trapped in co-rotation with the star by the stellar magnetic field. 
The clouds were first discovered as transient absorption features in Ho (Robinson t Collier 
Cameron 1986; Collier Cameron & Robinson 1989a: Paper I) . The absorption transients are 
produced when clouds pass between the observer and the stellar disc, scattering chromospheric 
Ha photons out of the line of sight. New clouds appear to form at the rate of one or two per 
day, just outside the Keplerian co-rotation radius, and subsequently move outVards until they 
finally escape from radial magnetic confinement, at least 10 R, from the rotation axis (Collier 
Cameron & Robinson 1989b: Paper II) . 

The specific angular momentum of the escaping cloud material is clearly high, and suggests 
that the cloud ejection process could contribute significantly to AB Dor’s rotational braking 
rate. Observations in Ha alone do not, however, provide very tight constraints on the cloud 
masses, since the population of the lower state of the transition depends critically on both 
the temperature and density in the cloud interior. The Call and Mgll resonance lines are 
much better suited to determining the cloud column densities, but time-resolved spectroscopic 
studies of the Call H &: K line profiles were not possible until dye-coated CCD detectors with 
extended blue response became widely available. 

In the present paper we report the results of a co-ordinated programme of observations of 
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AB Doradus made during 1988 December. We obtained high-dispersion, time-resolved Ha and 
Call H & K spectra at the European Southern Observatory (ESO) and the Anglo- Australian 
Observatory, and high-dispersion Mgll h & k spectra with the International Ultraviolet Ex- 
plorer ( IUE) satellite. In Section 2 we summarise the observations and data reduction tech- 
niques used, and in Sections 3 and 4 we describe the data analysis methods and absorption 
line profile synthesis used to constrain the projected areas, column densities, temperatures 
and internal velocity dispersions of the clouds. The implications for the masses and internal 
structure of the clouds are discussed in Section 5. 

2 Observations and data reduction 

2.1 Observing strategy 

The programme of IUE, A AT and ESO observations of AB Dor undertaken during 1988 De- 
cember is summarised in the journal of observations (Table 1). Of the four nights allocated to 
the programme at the AAT (Dec 16 to 19) only one (Dec 18) was lost due to poor weather, 
although the data quality during first half of the night of Dec 17 w r as degraded by variable cloud 
cover. The two ESA IUE shifts on 1988 Dec 17 and 19 coincided with the hours of darkness 
at the AAT Within the limits imposed by spacecraft temperature constraints and camera 
readout and preparation, the high-resolution LWP exposures of the Mgll h & Jt region were 
made at stellar rotation phases where well-defined absorption transients or clear patches had 
been observed in Ha at the AAT on Dec 16 and 17. Call H & K observations were interleaved 
with the Ha observations on Dec 19, the last night of the AAT run. 

All three nights of the ESO programme (Dec 21 to 23) were clear. The first night was 
devoted to continuous monitoring of the Ha profile, and the second to continuous monitoring 
of the Call H & if region. On the last night of the ESO run, Ha and Call H & K spectra were 
interleaved. During the three US IUE shifts on Dec 24 and 25, alternating LWP and SWP 
exposures were made, covering the full stellar rotation cycle. 
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2.2 A AT instrument configuration 

The A A T spectra were obtained with the new UCL echelle spectrograph (UCLES: Walker k 
Diego 1985) , at the coude focus of the 3.9m Anglo- Australian Telescope. The 31 lines mm -1 
echelle grating was used with the 70cm camera and a liquid nitrogen-cooled, dye-coated GEC 
CCD as detector, giving a dispersion of 2.7km s -1 per 22 fim pixel in the dispersion direction. 
A slit width of 1.3 arc sec gave an instrumental profile with full width at half maximum 
intensity of 2.8 pixel, or 7.5 km s During times of poor seeing, a focal reducer (Diego k 
Walker 1985) was used to increase throughput by opening the slit width to 1.7 arcsec, while 
maintaining the same resolution at the detector. The CCD chip was binned by a factor of four 
in the slit direction, to reduce both readout noise and chip readout time. The exposure times 
used were 150s for all Ha exposures, and 300s for Call H k K. 

2.3 ESO instrument configuration 

The ESO spectra were obtained with the Cassegrain Echelle Spectrograph (CASPEC; Pasquini 
& D Odorico 1989) on the 3.6m telescope at La Silla. A 31.6 lines mm -1 grating was used on 
the first night (1988 Dec 21 UT) in the red range centered at 6250.4. A 52 lines mm” 1 grating 
was used on the second night in the blue range centered at 420017 On the third night the 
31.6 lines mm -1 was used again, alternating between the blue and red ranges. The CASPEC 
echelle spectral range covers around 1000 A. at once, allowing features other than the prime 
lines (Ha and Call H k K) to be observed. The light is focussed by a 291mm camera onto a 
Thomson CCD detector with 1024 x 660 15 micron pixels. A slit width equivalent to 1.5 arcsec 
allowed an instrumental spectral resolution of 30000 or 10 km/s. The seeing ranged between 1 
to 1.5 arcsec and the weather was photometric during the ESO observing period. The whole 
CCD chip was read unbinned, using an exposure time of 300s for all AB Dor frames, with a 
time lag of 50-60s between exposures. The pixel to pixel response and wavelength calibration 
were determined by using a quartz flat-field and a Th-Ar lamp. 

2.4 IUE instrument configuration 

IUE (Boggess et al. 1978) high resolution spectra were taken through the large aperture (5 x 
10 arcsec) with the Long Wavelength Prime camera covering the Mgll h k k emission lines. 
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2.5 Data reduction: Ho and Ca II spectra 
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Data reduction: Mg II spectra 
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Here F C (X) is the continuum flux, which was fit as a second order polynomial. A c is the 
amplitude of the emission gaussian. F(A a ) is the value of the underlying stellar monochromatic 
flux at the wavelength A a of the interstellar absorption line, and (ignoring the contribution of 
the continuum) is given approximately by F(X a ) = A e exp[-z|(A a )]. W a is the equivalent 
width in angstroms of the interstellar line absorption and <j a is the width of the instrumental 
profile, treated here as a gaussian. We also use the notation 


z c (A) 

*-(A) 


A - A e 
A - A a 


( 2 ) 

( 3 ) 


where A c is the wavelength of the emission gaussian. As the equivalent width of the interstellar 
absorption and the instrumental resolution must be the same for all exposures, the lines from 
all exposures should be fit simultaneously. In practice, all of the k lines were fit simultaneously, 
and all the k lines were fit simultaneously, but the whole set of lines were never fit as one batch. 
The equivalent widths of the interstellar lines and the instrumental resolutions from the two 
different fits were quite similar. 

The fitting was done w T ith a modified amoeba algorithm (Press et al. 1986; Kuntz 1990, 
in preparation). Uncertainties determined for the calculated quantities are the most extreme 
values of those quantities for which the chi square of the fit (at the extreme value of the 
calculated quantity) minus the best possible chi square for the fit is unity. The algorithm used 
to determine uncertainties is described in full in Kuntz (1990, in preparation). 

The fluxes from the Mgll lines were measured in two ways: (a) the area under the fitted 
gaussian emission profile was calculated, and (b) the absorbed flux F(X a )W a calculated from 
the fit was added to the observed flux minus the fitted continuum, integrated over a 3. 4 A (i.e. 
±180 km s l ) window centered on the line. The measured flux (minus fitted continuum), the 
fitted absorbed flux, the corrected measured flux and the fitted emission flux are listed for each 
spectrum in Table 2. The two methods do not agree as well as one would like. The ratio of 
the quantity as measured by the first method to that measured by the second will always be 
greater than unity, and will increase with increasing noise. This is due to the fitted flux being 
a function of the width of the emission gaussian, which will always be larger than it should be 
due to the noise in the wings. The quantity derived by the second method (labelled “Corrected 
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flux” in Table 2) is therefore to be preferred. 


3 Measuring the cloud absorption profiles 

3.1 Ha profile analysis 

Normally, chromospheric emission fills the Ha profile in AB Dor, almost exactly to the level 
of the adjacent continuum. To a good approximation we can thus treat the background illu- 
minating source for the clouds as a flat continuum in Ha. When a cloud passes between the 
stellar disc and the observer, it scatters chromospheric Ha photons out of the line of sight. 
There is little coupling between the thermal properties of the cloud and the radiation field. In 
Paper I we inferred from statistical equilibrium studies that the cloud source function consists 
solely of the stellar radiation field, diluted according to the distance of the cloud from the 
star. The population of the n = 3 level of HI is sufficiently low that both recombination and 
stimulated emission in Ha are negligible. Since most of the clouds are located approximately 
3 R* from the star, the cloud source function should be at most 3 percent of the intensity of 
the chromospheric background in Ha. For this reason we expect to see a "pedestal” of weak 
emission coining from off-disc clouds. In Paper I we found that the amplitude of the pedestal 
emission in a time-averaged spectrum was greatest at ±3v e sin i from line centre, and was at 
most one or two percent of the stellar background spectrum in HaThis agrees with our expec- 
tations, because the majority of the clouds are located just outside the co-rotation radius and 
the clouds near the inner edge of the distribution will have the greatest surface brightnesses. 
The upper limit of 1 or 2 percent on the pedestal emission tells us that the total projected 
area of all the clouds in the system is probably between 1 and 5 times the projected surface 
area of the star. Since the average cloud cross-sectioned area is roughly 20 percent of the disc 
area, this suggests that at any time there are between 5 and 20 clouds of this size surrounding 
the stax, and that they form preferentially close to the stellar equatorial plane. 

The radial temperature gradient in the layers of the stellar chromosphere where Ha and 
the Call and Mgll resonance line emission features are formed is small enough that we are 
justified in neglecting centre-to-limb variations in the emission line profiles for a lowest-order 
analysis. 
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In this very simple situation the observed equivalent width IV of the cloud absorption 
profile is just the average of the absorption equivalent width w calculated along all lines of 
sight intersecting the stellar disc. Since the clouds appear to be well-defined, sharp-edged 
entities, we can also express the observed equivalent width in terms of the effective fraction of 
the stellar disc covered by clouds, and an area-weighted mean absorption equivalent width w: 


W 


<hi! Lc w(x ' y)iyix 

jij '^cloud 
A, ' 


( 4 ) 

( 5 ) 


If we can determine the value of w from the line profiles, we can use the observed EW variations 
in Ha to determine the fraction of the stellar disc covered by clouds at any time. 

In interpreting the shapes of the line profiles, we make the simplifying assumptions that 


1. the cloud source function is zero, 

2. the background continuum intensity T 0 is independent of wavelength, 

3. the absorption equivalent width w is the same along all lines of sight intersecting a cloud 
and the stellar disc, and 


4. there is no centre- to- limb variation in the chromospheric Ha surface brightness. 


When no cloud is obscuring the disc, the observed continuum flux is 


F c = Jq J c{x)dx = irl 0 


( 6 ) 


Here we define the independent variable x as the radial velocity along a given line of sight in 
units of v e sin i, or equivalently as the projected x-coordinate of the line of sight on the stellar 
disc in units of the stellar radius. A chord across the stellar disc, running through x parallel to 
the projected rotation axis, has a length c(x) = 2(1 - x 2 ) 1 / 2 . We also define a function h(x) to 
describe the shape of a section of the star obscured by cloud. The units of h(x) are the same 
as those of c(x), so that 

f h(x)dx = w — £ l ° u d . ( 7 ) 

J —\ Am 

The rotation profile of the star is thus given by c(x), and the rotation profile of the part of the 
stellar disc occupied by clouds is given by h(x). 
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The observed line profile is the convolution of the background spectrum I(x) = Jo with 
the rotation profile of the unobscured part of the disc, plus the convolution of the absorption 
profile J c (z) along a single line^of sight through the cloud with the cloud rotation profile: 

F(x) = Jo J j [c(x') - A(z')] dx' + J J c { x - x')h(x') dx'. (8) 

We obtain the fractional depression of the observed monochromatic flux below the level of the 
continuum from Eqs. (6) and (8): 

- = tt q ~ T <( x ~ X 'M X ') dx ' ( 9 ) 


In the limit where the cloud has a small extent in x y giving a rotation profile h(x) narrower 
than the width of the intrinsic absorption profile J c (x), we can treat h(x) as a ^-function with 
area 7Ti4 c j OU( j/A 1 ,, located at the x-coordinate x c of the cloud centre, Eq. (9) then simplifies 
to 

"" Fx _ ^cloud f jo ~ ~ J c) 

A* v Iq 

If the cloud absorption profile is saturated, then J c (0) = 0. The fractional depression of the 
deepest part of the observed line profile is then equal to the fraction of the stellar disc obscured 
by the cloud. Note that if we integrate over x } we recover the result of Eq, (5), since 



and 





Iq 


(ii) 


( 12 ) 


In practice, this condition is satisfied for a time as the leading edge of a cloud starts to 
move on to the disc, and again as its trailing edge moves off the disc. During these opening and 
closing stages of a transient, the equivalent width of the observed profile is W — ^^cioud/^*’ 
and the fractional depth of the centre of the absorption profile is (1 — F(x c )/F c ) = A^^ /A*. 
That is, a plot of the equivalent width against the depth of the observed absorption profile 
should give a straight line with slope w. 

Fig. 2 shows the width-depth plots for the opening and/or closing phases of twelve of the 
best-defined Ha absorption transients observed during the campaign. The slopes of the linear 
parts of these relations are listed in Table 3. The average value is 44.98 pixel (2.66A) with a 
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2a scatter of 12.56 pixel (0.74A). The small observed spread in the values of w allows us to 
estimate reliably the rate at which the fraction of the disc obscured by clouds is varying, even 
if we cannot make a direct density estimate for each individual cloud. 

3.2 Ca II profile analysis 

At the temperatures and densities inferred for the clouds in Paper I, most of the calcium is 
expected to be singly-ionised and in the ground state. In this case we can estimate the total 
column density of the clouds from the column density of Ca II ions, assuming normal solar 
abundances. 

The cores of the Call H k K resonance lines in AB Dor show strong chromospheric emis- 
sion reversals, giving a plentiful background source of photons in these lines. The effect of the 
circumstellar clouds on the observed profiles is similar to that in Ha: Call H k K photons are 
scattered out of the line of sight when clouds pass between the stellar disc and the observer, 
giving rise to travelling absorption features in the dynamic spectra. No strong off-disc emis- 
sion features are present in the dynamic spectra, indicating that once again the cloud source 
function in these lines is small, consisting mainly of scattered stellar Call H k K photons. 

In interpreting the shapes of the Call H k K profiles, we make simplifying assumptions 
similar to those used in the discussion of the Ha profiles: 

1. the cloud source function is zero, 

2. the total cloud column density is the same along all lines of sight intersecting a cloud 
and the stellar disc, and 

3. there is no centre-to-limb variation in the chromospheric Call H k K spectrum. 

If the underlying spectrum in the co-rotating frame is 7o(x), then the emergent spectrum along 
a line of sight through the cloud is 


Ic(x) = Io{x)e T W. (13) 

The observed monochromatic flux is obtained by convolving the underlying spectrum with 
the rotation profile of the unobscured part of the disc, and the absorbed spectrum with the 
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rotation profile of the obscured part: 


F(x) = J J 0 (x - x , )[c(x / ) - h(x')] dx* + J I c (x - x')h(x') dx f , (14) 

The total flux scattered out of the line of sight by the clouds is most conveniently expressed 
as a line index Ik or J//, measured relative to the total (unobscured) flux over a wavelength 
band extending from 1.5u c sin i blueward to 1.5u c sin i redward of line centre: 

J%[Fb(x l )-F(x')]dx' 


Ik.H = 


1% Fo(x')dx> 


(15) 


If we define a similar line index for a single line of sight through the cloud 

X- 3/2 [M*') - T c( x ')] dx‘ 


l K,H = 


f?i% Io{x’)dx' 


(16) 


we can see that to a good approximation, Ik,h = the case with the Ha 

equivalent widths. A plot of the observed line index Ik or Ik against Ha equivalent width W 
should give a straight line with slope i/w. This gives us a measure of 7, since we can determine 
w directly from the Ha profile width/depth variations: 


ik,h - xv 


dlK,H 

dW 


(17) 


In practice, this technique works best when we use short sequences of interleaved Ha 
and Call H h K observations made at times when the total amount of cloud coverage on the 
stellar disc is changing rapidly. Over longer periods, rotational changes in the background 
chromospheric emission flux can give rise to spurious drifts in the Ha equivalent width and 
the Call H ic K line indices. 

For shorter sequences, the dominant source of error in measuring the line core fluxes in 
the calibrated, normalised spectra lies in the polynomial fit used to straighten the continuum. 
We compensate for this by measuring the line core flux relative to two reference bands, each 
0.75u e sin i wide, immediately adjacent to the line core band. Photon shot noise and readout 
noise then provide a good estimate of the uncertainty in each line index measurement. For the 
Dec 19 and 23 data, we interpolate linearly in time to obtain the value of W at the time of 
each Ca n observation. For the continuous series of Ca II spectra taken on 1988 Dec 22, we 
can estimate the the values of W at a given rotation phase by interpolation, using Ha spectra 
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taken at the same rotation phase on 1988 Dec 21 and 1988 Dec 23. Fig. 3 shows the difference 
between values of W measured on these two nights, at the phase corresponding to a given 
UT on 1988 Dec 22. During the intervals used in Fig. 4 and Table 5, it can be seen that the 
phase*dependent changes in W did not evolve significantly between Dec 21 and Dec 23. Note, 
however, that the time-dependent zero-point drifts in the Ik,H - W relations for Dec 22 are 
considerably greater than those on Dec 19 and 23. The Ha equivalent widths obtained using 
this technique may not be as reliable as those obtained by alternating observations at the two 
different wavelengths. 

In Fig. 4 we show the Call H Sz K line index - Ha equivalent width relation for a total of 
8 sequences on Dec 19, 22 and 23. The slopes of the straight-line fits to these plots are listed 
in Table 5. 

3.3 Mg II profile analysis 

The analysis of the Mgll h & k profile variations on 1988 Dec 17 and 19 is very similar to that 
described above for the Call H & K data. Like the calcium, most of the magnesium in the 
clouds will be singly-ionised and in the ground state. The chromospheric Mgll h & k emission 
provides a background source of photons which are scattered out of the line of sight by the 
clouds. Again the line source function in the clouds is expected to consist solely of scattering 
terms. The low level of the surrounding continuum allows us to place further constraints on the 
total projected surface area of the cloud system. The base widths of the Mgll h & k emission 
lines are close to 3Aclouds located just outside the co- rotation radius should produce pedestal 
emission between 2.0 and 2.5A to either side of line centre when they are near maximum 
elongation. The pedestal emission comprises no more than 30 percent of the signal at this 
distance from line centre, and may be less because of the extended wings of the chromospheric 
emission features. If, as we argue below, the cloud optical thickness is such that the width 
of the absorption profile is similar to the emission base width, then the total scattered flux 
from the whole cloud system should be some 3 percent of the total chromospheric emission 
flux multiplied by the total projected area of the cloud system in units of the stellar disc area. 
The upper limit on the total scattered flux is of order 4xl0~ 13 ejgcm 2 s* 1 , while the total 
chromospheric emission flux is of order 4x 10~ 12 erg cm 2 s” 1 . This implies a total projected 
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area for the cloud system of three times the stellar disc area, in good agreement with the value 
inferred from the Ha pedestal emission. 

The procedures for forming line indices and calculating their dependence on cloud column 
density are nearly identical to those described above for Ca II, and so will not be repeated 
here. The band over which the line indices I k j t and , are calculated extends from 2.0 v e sin i 
blueward to 2.0 v c sin i redward of line centre. No photospheric reference sidebands are used 
in this case because of the very low level of the surrounding continuum. The sensitivity of 
the LW P camera is sufficiently stable that the use of such sidebands is unnecessary. The 
line profiles and integrated fluxes have, however, been corrected for the effects of interstellar 
absorption as described in Section 2.6 above. 

The main problem in interpreting the simultaneous Ha and Mgll h k k data is the long ex- 
posure times of the IUE spectra. Substantial changes in cloud covering fraction often occurred 
during the Mgll exposures, as seen from the changes in the Ha EW during each exposure. A 
second problem is that it was not possible to measure the intrinsic Ha EW (uf) for every cloud 
contributing to the Mgll line index - Ha absorption EW relations. 

The best we can do is to average all measurements of the Ha absorption equivalent width 
made during each Mgll exposure, to obtain the time-averaged value of W during the IUE 
exposure. Secondly we must assume an average value of W in order to estimate the average 
cloud covering fraction durung the exposure from the average value of W. We must bear in 
mind that the results may be affected by spurious zero-point drifts, caused by rotationally- 
dependent changes in either the underlying Ha or Mg II emission flux. However, the Ha 
width-depth plots for five different transients observed on 1988 Dec 17 and 19 all lie on the 
same relation (Fig. 2). This suggests that any such drifts are small, in Ha at least. 

The observed Mgll h & k line indices and the corresponding time-averaged W values 
for the Ha absorption are listed in Table 6. The I kfh - Ha absorption EW relations are shown 
in Fig. 5. The discrepant point in both these plots corresponds to the first IUE exposure made 
on 1988 December 19, LWP14672. Closer examination of this spectrum revealed that both 
the Mgll emission and the surrounding continuum fluxes are consistently lower than the other 
exposures, by nearly a factor of 2. The reason for this flux deficiency is not clear: it is possible 
that the star was not properly centred in the entrance aperture. 
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If the data from LWP 14672 are excluded, the 7^ - Ha absorption EW relations are 
well-represented by linear relations, with slopes of (dI k /d\V) = 0.58 ± 0.06 and ( d.I h (dW ) = 

0.49 i 0.06 respectively. If we assume that the cloud scatters all the chromospheric Mgll 
photons out of any line of sight intersecting the cloud, this implies an upper limit on w of 
1.90 ± 0.15 A . This is within the scatter of the distribution of w values measured from the Ha 
width-depth relations. This indicates that the width of the absorption profile of a typical cloud 
in the Mgll resonance lines is comparable with the base width of the chromospheric emission 
features, and so places a lower limit on the cloud optical thickness. 

4 Line index - column density relations 

4,1 Ha profile synthesis 

For the larger clouds studied here, the range in intrinsic Ha absorption equivalent widths UJ 1 
and hence in column densities of hydrogen atoms in the first excited state (n = 2), is small 
(Table 3). These larger clouds have projected surface areas between 13 and 18 percent of the 
stellar disc area, so it is not too surprising that their other physical properties should be fairly 
similar. 

The equivalent width of the Ha absorption profile depends principally on 

1. the total column density of cloud material; 

2. the degree of excitation and ionisation of hydrogen in the cloud interior; and 

3. the amplitude of turbulent random motions within the cloud. 

turbulent motions within the cloud have the effect of broadening the Doppler core of the 
absorption profile, although it is less important for hydrogen than for heavier atoms with 
small thermal doppler widths. By analogy with quiescent solar prominences, we suspect that 
turbulent motions may be important in the AB Dor clouds. The clouds are likely to have a 
highly-fragmented structure: the effective pressure scale height in the cool cloud material is 
determined by the combined gravitational acceleration and the cloud temperature, and is of 
order 200 to 300 km, which is very much smaller than the overall cloud extent. Nakagawa 
iz Malville (1969) found that in the solar case, the prominence-corona interface was unstable 
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to fragmentation for perturbations with wavelengths longer than the pressure scale height. 
Given that conditions inside the AB Dor clouds are likely to be similar, it seems plausible that 
filamentary structures similar to those seen in quiescent solar prominences may be present, 
and that individual elements of this fine structure may have random velocities at least as great 
as the 5 to 10 kms 1 observed by Engvold (1972) in the solar case. 

In Fig. 6 we show contours of constant Ha absorption equivalent width W, computed as a 
function of the column density of H atoms in the n = 2 level and a gaussian turbulent velocity 
component. The contours span the observed range of w values, and show little change as the 
turbulent velocity dispersion is increased from 0 to 20 kms -1 . The range in column density of 
singly-excited H atoms permitted by the observations spans less than an order of magnitude: 
for a turbulent velocity of 10 kms -1 , logn n=2 = 14.4 ± 0.2cm -2 . 

The degree of excitation of hydrogen within the cloud is influenced mainly by the cloud 
kinetic temperature, and to a lesser extent by the volume density of the cloud material. The 
statistical equilibrium in the cloud interior was discussed in some detail in Paper I, to which 
the reader is referred for further details. Here we use the same model with a 5-level H atom 
plus continuum, and the Lyman lines and continuum in radiative detailed balance. In Fig. 7 
we show contours of constant in as a function of cloud temperature and total hydrogen column 
density n// + ftp, computed using a turbulent gaussian velocity dispersion of 10 kms -1 . 

Fig. 7 was computed for a slab of uniform density with a thickness of one stellar radius, 
comparable to the projected linear dimensions of the clouds. If the clouds are significantly 
fragmented, however, the local density in the neutral material will be very much greater than 
we have assumed. We can explore this possibility by computing a second set of models in which 
the slab thickness is decreased by a factor of 1000 and the local density increased by the same 
factor to give the same total column density. The results axe shown in Fig. 8. At temperatures 
below 8000 K there is little difference between the uniform and fragmented models. At higher 
temperatures, however, when a substantial fraction of the hydrogen is ionized, the higher 
local densities in the fragmented model serve to increase the recombination rate, increasing 
the amount of neutral hydrogen relative to the uniform density model. The observed optical 
thickness in Ha is thus achieved at a column density up to an order of magnitude lower in the 
fragmented model than in the uniform model. 
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Figs, 7 and 8 show clearly that the observed Ha equivalent widths effectively rule out total 
cloud column densities less than (8 ± 4)xl0 18 cm“ 2 . This lower limit on the column density 
occurs at a temperature close to l.lxl0 4 K, where the number density ratio of H atoms in 
the n = 2 level passes through a maximum relative to the total number density (n// + np). 
At higher temperatures, the increasing degree of excitation is offset by the decrease in the 
amount of neutral hydrogen as the hydrogen becomes ionised; while at lower temperatures the 
hydrogen is mainly neutral but the degree of excitation is lower, 

4.2 Ca II profile synthesis 

We now need to calibrate the variation of ik,H with column density. The amount of flux 
scattered out of a single line of sight depends on ^ 

1. the width of the chromospheric emission core in the background stellar spectrum as seen 
in the co-rotating frame of the cloud ; 

2. the velocity dispersion of the Ca II ions in the cloud; and 

3. the column density of Ca II ions in the cloud. 

Although we cannot determine the intrinsic width of the emission core directly (owing 
to the extreme rotational broadening of the spectrum) we can make a reasonable estimate 
based on observations of more slowly-rotating, active stars of similar spectral type. Empirical 
studies of Call H & K emission profile widths in a large number of dwarfs and giants with 
a range of spectral types and chromospheric activity levels show that the full widths at half 
maximum (FWHM) of the chromospheric Ca II emission cores are insensitive to the level of 
chromospheric heating, although they are very sensitive to surface gravity (Wilson & Bappu 
1957) . In a theoretical study of chromospheric scaling laws, Ayres (1979) found that the base 
width A A* of the chromospheric Ca II emission cores and the thermalisation width AA^ have 
equal and opposite dependences on met alii city, stellar effective temperature and chromospheric 
heating, but the same dependence on surface gravity. The thermalisation width corresponds 
to the separation of the red and blue peaks of the emission profile. The functional dependence 
of the FWHM of the emission profile on the various stellar parameters is intermediate between 
those of the base width and the thermalisation width: that is, the effects of metallidty, effective 
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temperature and chromospheric heating cancel out, leaving only the gravity dependence of the 
Wilson-Bappu effect. 

The compilation of high-resolution Call H K profiles of late-type stars by Linsky et al. (1979) 
gives FWHM values for several stars of similar surface gravity and spectral type to AB Dor. The 
scatter in these widths is small, and w r e adopt a mean value of 0.43±0.04A as the FWHM of the 
Call H & K emission cores in AB Dor. For modelling purposes, we represent the background 
spectrum by a gaussian emission core of appropriate width and amplitude, superimposed on a 
parabolic fit to the broad photospheric absorption profile. 

The atomic parameters used in calculating the line absorption profiles are listed in Table 4. 

The thermal Doppler width of the absorption profile is 1.2±0.2kms“ 1 over the temperature 
range 5000 to 10000 K, in which the cloud temperature is expected to lie. As described above, 
however, we also consider the possibility that turbulent broadening of up to 20 km s'" 1 may be 
also be present due to random motions of fragmented sub-structures within the clouds. 

In Fig. 9 we plot contours of constant line index ix and in as a function of cloud column 
density and turbulent velocity, for an assumed cloud temperature of 8500 K and a background 
emission core width of 0.43A. The effects of varying the background emission core width are 
illustrated in Fig. 10: the effects on the line index are negligible. 

We can use the model line index maps shown in Fig. 9 to determine a confidence region 
in (n// + rip, u tur b) space, within which we can achieve a satisfactory fit to the observations. 

For each sequence of n Ca II observations of a single absorption transient, we usually had 
observations of the K-line in two different echelle orders, and of the H-line in a single order. 

Each ( nx + n p , r tur b) pair defines the slope of a theoretical line index - Ha equivalent width 
relation, in each of the H and K lines. At each value of n// + n p and u tU rb fix the slope 
at the theoretical value, and determine the zero-point offset for each line in each order, by 
minimising the statistic 


5 = Yl n 


Di - Fi 


1 = 1 


Vi 


(18) 


where Z), and Fi are the observed and predicted data values and < 7 , is the observational scatter. 
We then sum the resulting values of 5 for all lines and orders, and so build up a map of 5 as 
a function of n// + n p and r tU rb> incorporating all the data on each individual cloud observed. 
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We identify the contours of constant 5 corresponding to the 95.4 and 99.0 percent confi- 
dence limits on (n// + n p , i> lur b) using the parameter estimation method of Lampton, Margon 
k Bowyer (1976). The value Si on the contour that has a probability a of enclosing the true 
values of the parameters is given by 

S L = Smin + X p (a). (19) 

Here 5min the minimum value of 5* in the map, and x p (^0 is the tabulated value of the 
X 2 distribution for p degrees of freedom and significance a. In this case there are two free 
parameters, so that p - 2. The 95.4 and 99,0 percent confidence regions in (n// + n p , u tur b) 
space are shown for the eight clouds observed, in Fig. 11. 

We verified that p = 2 was the correct choice by Monte Carlo simulation, as follows. We 
assume that the “true” cloud parameters are those giving the minimum value of 5, i.e. 5 truc = 
‘S'min- Using these parameters we generate a large number of artificial data sets with normally- 
distributed errors. For each artificial data set we generate a new map of 5 / (n// + n p , Vturb)> and 
locate the minimum value S^ n . The value of 5 corresponding to this (n H + n p , v turb ) pair in 
the original map derived from the actual data is noted. After n such simulations, we sort the 
n resulting S values into ascending order, and note the limiting values Si below which 95.4 
and 99,0 percent of the simulated minima lie. We found good agreement between the contour 
levels estimated in this way and those estimated by the method of Lampton et al 

4,3 Mg II profile synthesis 

The modelling of the Mg II line profiles as viewed in the co-rotating frame (and hence the Mg 
II line indices ih t k) as functions of the cloud column density was very similar to that used for 
the Ca II profiles. Again the chromospheric emission profile was approximated by a gaussian. 
The width of the synthetic profile was based on IUE observations of chromospheric Mg II 
emission profile widths for several nearby late-G and early-K dwarfs (Basri k Linsky 1979). 
As with the Ca II profiles of Linsky et al. (1979), the scatter in the profile widths for the stars 
in the Basri k Linsky sample is small, and we adopt a mean FWHM of 0 . 60 ± 0 . 06 A for our 
synthetic gaussian profile. 

Fig. 12 shows contours of constant line index i * and ih as functions of cloud column density 
and turbulent velocity, for an assumed cloud temperature of 8500 K and a background emission 
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core width of 0.60A. In Fig. 13, the effects on the line index - column density relations of varying 
the background emission core width by ±0.06Aare represented as in Fig. 10: once again, the 
line index is seen to be insensitive to errors in the estimated width of the chromospheric 
emission core. 

The contours of constant x 2 bounding the 95.4 and 99.0 percent confidence regions in 
(n// + n p , t’turb) space, within which we can obtain a satisfactory fit to the combined 1988 
Dec 17 and 19 Mgll h &: k line index - Ha EW relations, excluding LWP 14672, are shown in 
Fig. 14. 

5 Cloud masses and internal velocity fields 

Using the techniques developed in Section 3 we have succeeded in measuring both the equivalent 
width of the Ha absorption and the amount of the Mg II and Ca II resonance line scattering, 
along typical lines of sight running through some of the larger circumstellar clouds. When 
compared with the line profile synthesis results of Section 4, the observed line indices provide 
a set of independent constraints on the column densities, temperatures and turbulent internal 
velocity dispersions of the clouds. The Ha equivalent widths are more sensitive to the density 
and temperature of the cloud material than to the random internal motions. The Call H fz K 
and Mgll h & k line indices are more sensitive to column density and turbulent motions than 
to cloud temperature. 

The Ha equivalent widths of the various clouds studied show remarkably little scatter. If 
the physical conditions affecting the degree of excitation and ionization of hydrogen are similar 
for all the clouds, this would imply that the total hydrogen column density is similar for the 
various clouds studied. The situation for the Ca II resonance lines is not quite so simple. 
For a fixed turbulent velocity dispersion, the scatter in the inferred column density varies by 
as much as two orders of magnitude from cloud to cloud. The three clouds for which w r e 
obtained good Ca II line index - Ha equivalent width relations on 1988 Dec 22 and 23 showed 
a marked decrease in Ca II absorption over the course of two stellar rotations. Over the same 
period, however, the Ha equivalent width and fractional area coverage underwent only minor 
evolution. 

There are two possible explanations for this behaviour. First, the amount of material in 
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each of these three physically distinct clouds may have decreased by a factor of between 10 and 
100 between 1988 Dec 22 and 23. Let us suppose that the clouds are in thermal equilibrium, 
such that a fixed rate of conductive and mechanical energy input is balanced by radiative 
losses in strong, effectively thin lines such as the hydrogen Balmer series and the Mg II and 
Ca II resonance lines. In this case, a drop in the cloud density will lead to an increase in 
temperature until the electron density (and hence the radiative loss rate) has attained more 
or less its original value. At the temperatures and densities considered here, the n — 2 level 
population of hydrogen varies almost in direct proportion to the electron density, provided that 
the degree of ionization of the cloud material is small. This “thermostat” acts to keep the Ha 
absorption equivalent width at a nearly constant value. At temperatures above 10 4 K, however, 
hydrogen is almost completely ionised so that this mechanism breaks down. As Figs. 7 and 8 
show, the lowest total column density that can possibly give rise to the observed Ha equivalent 
widths is of order 10 39 cm* 2 , so that the total column densities on 1988 Dec 22 must have been 
one to two orders of magnitude greater, i.e, at least 10 2 ° to 10 22 cm" 2 . 

The second possibility is that the column density and temperature did not change. Instead, 
the decrease in the strength of the Ca II absorption could have been caused by a change in 
the average amplitude of the turbulent motions inside the clouds, from 10-20 km s“ 3 on Dec 
22 to 4kms" 3 or less on Dec 23. The solar analogy provides some evidence in favour of this 
idea. Solar quiescent prominences often undergo temporary increases in the amplitude of their 
internal motions. These activations usually occur in response to a solar flare, when large-scale 
magnetic wave disturbances spread out through the corona, propagating large distances over 
the surface of the Sun at speeds of order 1000 km s" 1 . The internal motions in the prominence 
increase from the usual 5-10 km s" 1 to 30-50 km s" 1 , and after an hour or so may either die 
down or lead to an outward eruption of the entire prominence (Bruzek 1969). Given that AB 
Dor undergoes strong soft X-ray flares at a rate of one or two per day (Pakull 1981; Vilhu k 
Linsky 1987; Collier Cameron et aL 1988), this explanation seems plausible and has the added 
attraction that it can explain why three physically distinct circumstellar clouds should all show 
enhanced Ca II absorption only during one four-hour series of observations, on Dec 22. Further 
simultaneous optical and X-ray observations are planned for the near future, and should allow 
us to determine whether X-ray flares and Ca II absorption enhancements are correlated. 
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In Table 7 we list the various constraints on the cloud temperature and internal motions, 
as functions of column density. The range of cloud temperatures consistent with the observed 
Ha equivalent width is listed, for a cloud model with uniform density (filling factor / = 1, as 
in Fig. 7) and for a fragmented model with filling factor / = 10" 3 (Fig. 8). The 99 percent 
confidence limit on the value of t; tur b consistent with the observed amount of Mg II absorption 
on 1988 Dec 17 (Fig. 14) is shown, as are the limits on ^turb for each of the clouds observed 
in Call H & K (Fig. 11). Under the assumption that the differences in Ca II absorption 
from cloud to cloud are due to differences in internal velocity dispersion rather than column 
density, it seems that the only value of the column density that is capable of satisfying all 
the observational constraints simultaneously, is nn + n p ~ 10 20 . The corresponding cloud 
temperature needed to give the observed Ha equivalent width lies in the range 8000 to 9000 K, 
and the internal motions of the clouds range fronTless than 2kms~ 1 to more than 14kms -1 . 

The cross-sectional areas of all the clouds studied in Call H & K and Ha were in the range 
0.15 to 0.20 times the projected area of the stellar disc. For a column density of 10 2O cm~ 2 , 
this gives a mass of 4 to 5xl0 17 g per cloud. The amount of energy needed to increase the 
turbulent velocity dispersion of such a cloud from ^5 to ^ISkms -1 is of order 5xl0 2 ^ erg - a 
small fraction of the ^10 34 erg released in 0.1 to lOkeV X-rays by a single large flare on AB 
Dor (Collier Cameron et al 1988).. 

The energy budget for disposing of between 90 and 99 percent of the material in three 
clouds between 1988 Dec 22 and 23 is considerably greater. A large fraction of the cool cloud 
material seen on 1988 Dec 22 material must either have been re-ionised, fallen back to the 
surface of the star or been ejected from the corona by Dec 23. The amount of energy needed to 
re-heat most of the material in a 4x 10 17 g cloud to the ambient coronal temperature of 2x 10 7 K 
is of order 10 33 erg. Alternatively, the confining field may have contracted, dragging the cloud 
down to a position inside the co-rotation radius, where the dense cloud material is free to drain 
back on to the chromosphere, releasing 10 32 to 10 33 erg of gravitational energy per cloud as 
it goes. A third possibility is that part of each cloud may have undergone centrifugal ejection 
from the confining field, again releasing at least 10 33 erg of magnetic energy from the stressed 
magnetic configuration. However, the apparent stability of the radial positions and projected 
areas of these three clouds, and the lack of any obvious bulk motions (absorption blueshifts or 
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redshifts) during the period 1988 Dec 21 to 23, tends to argue against any hypothesis requiring 
all three clouds to shed most of their contents in such a spectacular fashion in the course of 24 
hours. 

6 Conclusions 

The new techniques developed in this paper for measuring the strength of the Ha, Call H & K 
and Mgll h &: k absorption along lines of sight intersecting a cool circumstellar cloud have 
been applied successfully to new optical and UV spectra of AB Dor. These measurements 
place important new constraints on the masses and internal velocity dispersions of the AB 
Dor cloud system. The average amount of Ha absorption along lines of sight through several 
of the larger clouds shows little variation, either from night to night or from cloud to cloud. 
The Call H K absorption, however, can change dramatically from night to night in a given 
cloud. These changes correspond to changes of 1 or 2 orders of magnitude in the cloud column 
density at a fixed internal velocity dispersion. A more plausible explanation for these changes 
in the Call H & K absorption is that the internal velocity dispersions of the clouds may vary 
in time over the range 5 to 20kms“ 1 - possibly increasing in response to stellar flares - in 
a manner analogous to the activation of solar prominences. In this case, a column density 
n H + ftp = 10 2O cm~ 2 is consistent with all the observations. The projected cross-sectional 
areas of these clouds are 15 to 20 percent of the projected stellar disc area, giving total masses 
per cloud of 2 to 6xl0 17 g. 

The excitation temperature for hydrogen needed to match the observed Ha EW is between 
8000 and 9000 K. If we take into account the combined local gravitational and centrifugal 
acceleration, the pressure scale height in the clouds is thus a few hundred km. By analogy 
with solar quiescent prominences, where similar physical conditions prevail, we suspect that 
the cloud material may be concentrated in filamentary structures with local densities two 
to three orders of magnitude greater than the average cloud density. The amplitude of the 
turbulent velocity dispersion inferred from the Call If & K absorption strength is similar to 
the amplitude of the random motions of filaments within quiescent solar prominences. 

At any given time there appear to be of order 6 to 10 large clouds in the observable slice of 
the corona, with allowance made for incomplete phase coverage. The actual number of clouds 
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is probably at least twice this number, giving a total mass for the cloud system of order 10 19 
g- 

Acknowledgments ACC thanks the director and staff of the Anglo- Australian Observatory, 
for support at the telescope and for the use of the VAX 11/780 computer at AAO Epping, 
where the early stages of the data reduction were carried out. BHF and PE thank the staff of 
the European Southern Observatory for their support during the La Silla observations. MVP 
thanks the Institute of Astronomy, Cambridge for hospitality during the period over which this 
work was carried out, and the staff at Vilspa for support during the ESA IUE observations. 
The later stages of the data analysis were performed on the STARLINK VAX 11/780 at the 
Royal Greenwich Observatory. During the course of this work ACC was supported in a SERC- 
funded postdoctoral fellowship at the University of Sussex. An EEC grant to BHF for research 
into stellar activity is also acknowledged. 


25 


References 


[1] Ayres, T.R., 1979. Astrophys. J., 228, 509. 

[2] Basri, G.S. & Linsky, J.L., 1979. Astrophys. J., 234, 1023. 

[3] Bruzek, A., 1969. In de Jager, C. Sc Svestka, Z., editors, Solar Flares and Space Research , 
page 61, North-Holland, Amsterdam. 

[4] Collier Cameron, A., Bedford, D.K., Rucinski, S.M., Vilhu, O., & White, N.E., 1988. 
Mon. Not. R. astr. Soc., 231, 131-147. 

[5] Collier Cameron, A. & Robinson, R.D., 1989a. Mon. Not. R. astr. Soc., 236, 57. 

[6] Collier Cameron, A. & Robinson, R.D., 1989b. Mon. Not. R. astr. Soc., 238, 657. 

[7] Diego, F. & Walker, D.D., 1985. Mon. Not. R. astr. Soc., 217, 347. 

[8] Eggen, O.J., 1973. Publ. astr. Soc. Pacif., 85, 289. 

[9] Eggen, O.J., 1975. Publ. astr. Soc. Pacif., 87, 37. 

[10] Eggen, O.J., 1983. Mon. Not. R. astr. Soc., 204, 377. 

[11] Eggen, O.J., 1983. Mon. Not. R. astr. Soc., 204, 391. 

[12] Engvold, O., 1972. Solar Phys., 23, 346. 

[13] Innis, J.L., Thompson, K., & Coates, D.W., 1986. Mon. Not. R. astr. Soc., 223, 183. 

[14] Lampton, M., Maxgon, B., & Bowyer, S., 1976. Astrophys. J., 208, 177. 

[15] Lites, B.W., Skumanich, A., Rees, D.E., Murphy, G.A., & Carlsson, M., 1987. Astro- 
pbys. J., 318, 930. 

[16] Nakagawa, Y. & Malville, J.M., 1969. Solar Phys., 9, 102. 

[17] Pakull, M.W., 1981. Astr. Astrophys., 104, 33. 

[18] Pasquini, L. & D Odorico, S., 1989. CASPEC: ESO Operating Manual 2, European 
Southern Observatory. 


26 



[19] Press, W.H., Flannery, B.P., Teukolsky, S.A., & Vetterling, W.T., 1986. Numerical 
Recipes: The Art of Scientific Computing. Cambridge University Press, Cambridge. 

[20] Robinson, R.D. Sc Collier Cameron, A., 1986. Proc. Astr. Soc. Australia, 6, 308. 

[21] Rucinski, S.M., 1982. Inf. Bull. var. Stars, 2203. 

[22] Rucinski, S.M., 1985. Mon. Not. R. astr. Soc., 215, 591-614. 

[23] Shine, R.A. Sc Linsky, J.L., 1974. Solar Phys., 39, 49. 

[24] Stauffer, J.R. & Hartmann, L.W., 1987. Astrophys. J., 318, 337. 

[25] Stauffer, J.R., Hartmann, L., Soderblom, D.R^, Sc Burnham, N., 1984. Astrophys. J., 280, 
202 . 

[26] van Leeuwen, F. Sc Alphenaar, P., 1983. In Rodono, M. Sc Byrne, P.B., editors, IAU 
Colloquium 71, Activity in Red Dwarf Stars, Reidel. 

[27] van Leeuwen, F., Alphenaar, P., Sc Meys, J.M.M., 1987. Astr. Astrophys. Suppl., 67, 483 
-506. 

[28] Vilhu, O. Sc Linsky, J.L., 1987. Publ. astr. Soc. Pacif., 99, 1071. 

[29] Vilhu, O., Gustafsson, B., Sc Edvardsson, B., 1987. Astrophys. J., 320 , 850. 

[30] Walker, D.D. Sc Diego, F., 1985. Mon. Not. R. astr. Soc., 217 , 355. 

[31] Wilson, O.C. Sc Bappu, M.K.V., 1957. Astrophys. J., 125 , 661. 


27 


Figure 1: A set of representative Ha and Call H and K spectra. 

Each series of spectra traces the development of a typical cloud absorption transient. Panel 
(a) shows the development of a transient in Ha on 1988 Dec 21. Panels (b) and (c) show 
smoothed reference spectra of the Call H K lines. These reference spectra were subtracted 
from the observed spectra to give the time-series of Call H ic K difference spectra in panels 
(d) and (e). These were obtained over the same range of rotation phases as (a), on 1988 Dec 
22, and are plotted on the same vertical scale as (b) and (c). 


Figure 2: Ha absorption equivalent width plotted against central absorption depth. 

A total of 12 width-depth relations for individual clouds are shown, measured from the dynamic 
spectra of (a) 1988 Dec 16, (b) Dec 17, (c) Dec 19, (d) Dec 21 and (e) Dec 23. The slope of 
each plot gives the average Ha absorption equivalent width w for a single line of sight through 
the cloud. 


Figure 3: Ha absorption EW plotted against UT on 1988 Dec 22. 

The variation of Ha EW with rotation phase on 1988 Dec 22 was assumed to be the mean of 
the values observed at the same rotation phase on Dec 21 and 23. The dot-dashed curve shows 
the difference between the EWs observed on Dec 21 and 23. The horizontal solid line segments 
denote the three phase ranges over which the Call H & K index - Ha EW relations for Dec 
22 were calculated. Although some changes in the cloud distribution clearly occurred between 
Dec 21 and Dec 23, the night-to-night changes over these phase ranges were small compared 
to the amplitude of the phase- dependent EW variation. 



Figure 4: Call H and K line index plotted against Ha absorption equivalent width. 

Data are given for each of the Call H & K lines in each echelle order observed, for a total of 
eight transients observed on (a) 1988 Dec 19, (b) Dec 22 and (c) Dec 23. The slope of each 
plot gives the value of za^/z/uT for a single line of sight through the cloud. 


Figure 5: Mgll h and k line index plotted against Ha absorption equivalent width. 

The line indices for the five IUE LWP spectra obtained on 1988 Dec 17 are shown as open 
circles. The slope of each plot gives the average value of ik,h/ tv for a single line of sight through 
the cloud material obscuring the disc at the time of observation. 


Figure 6: Intrinsic Ha absorption EW as a function of n n=r2 and v tU rb* 



Figure 7: Intrinsic Ha absorption EW as a function of nu + n p and T c i ou d. 

The cloud thickness is 7.3xl0 10 cm and the turbulent velocity dispersion is u tur b = lOkms -1 . 


Figure 8: Intrinsic Ha absorption EW as a function of n# + n p and Td ou d. 

The cloud thickness is 7.3 xlO 7 cm and the turbulent velocity dispersion is u tur b = lOkms -1 . 
This simulates a fragmented cloud, in which the condensed cloud material occupies only 10“ 3 
times the total cloud volume. 


Figure 9: Call H and K absorption line index map. 

Contours of constant (a) iff and (b) iff are plotted as a function of toted hydrogen column 
density tih + n p and turbulent velocity dispersion u tur b. From bottom to top, contour levels 
are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. 


Figure 10: Call H and K line index plotted against total hydrogen column density 
The solid curve shows the variation of the line index with column density, for zero turbulent 
velocity dispersion and a gaussian chromospheric emission profiles with FWHM = 0.43A. The 
effects of increasing or decreasing the chromospheric emission width by 0.04A are shown by the 
dashed and dotted curves respectively. The effect of introducing a turbulent gaussian velocity 
dispersion of 4kms _1 is shown by the dot-dashed curve. 


Figure 11: Confidence limits on cloud parameters from Call H and K data. 

The 95.4% and 99% confidence limits on the cloud parameters are shown in the (riff + n p , 
^turb) plane. The confidence limits are contours of constant x 2 , for all the data available on 
each of the eight individual absorption transients studied. 


Figure 12: Mgll h and k absorption line index map. 

Contours of constant (a) ** and (b) ih are plotted as a function of total hydrogen column 
density rc// + n p and turbulent^ velocity dispersion r tU rb- From bottom to top, contour levels 
are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, CL7, 0.8, 0.9. 


Figure 13: Mgll h and k line index plotted against total hydrogen column density. 

The solid curve shows the variation of the line index with column density, for zero turbulent 
velocity dispersion and a gaussian chromospheric emission profiles with FWHM = 0.66A. The 
effects of increasing or decreasing the chromospheric emission width by 0.06A are shown by the 
dashed and dotted curves respectively. The effect of introducing a turbulent gaussian velocity 
dispersion of 4kms’ 1 is shown by the dot-dashed curve. 


Figure 14: Confidence limits on cloud parameters from Mgll h and k data. 

The 95.4% and 99% confidence limits on the cloud parameters are shown in the (n// + n p , 
rWb) plane. The confidence limits are contours of constant x 2 , for a fit to the Mgll h iz k 
data obtained on 1988 Dec 17. 


Table 1: Journal of observations 


UT Date 

Ground-based optical 

IUE 

1988 

Ha 

Call Hk K 

Mgll h k k 

Dec 16 

AAT 12:22 - 17:40 

— 



Dec 17 

AAT 10:15 - 17.56 

— 

ESA 09:51 - 16:40 

Dec 18 

AAT 16:29 - 17:56 

— 

— 

Dec 19 

AAT 10.42 - 18:04 

AAT 11:58 - 17:42 

ESA 13:19 - 16:17 

Dec 21 

ESO 01:29 - 08:40 

— 

— 

Dec 22 

— 

ESO 01:50 - 07:45 

— 

Dec 23 

ESO 02:49 - 08:01 

ESO 03:08 - 08:07 

— 

Dec 24 

— 

— 

— 

Dec 25 

— 

— 

US2 01:16 - 08:11 

Dec 25 

— 

— 

US1 17:17- 23:11 

Dec 26 

— 

— 

US2 00:29 - 05:59 


Table 2: Mgll h and k emission flux measurements. All fluxes are in units of 10“ 12 ere cm 2 s -1 . 

TTITTk T ' •» r < . . . _ . ° 


LWP No 

Line 

Measured 

Absorbed 

Corrected 

Fitted 

Ratio 



Flux 

Flux 

Flux 

Flux 


14655 

k 

2.97±0.14 

0.92±0.02 

3.88±0.14 

4.87±0.12 

1.25 

14655 

h 

3.46±0.09 

0.75±0.01 

4.21±0.10 

4.66±0.09 

1.10 

14656 

k 

2.18±0.13 

0.63±0.02 

2.82±0.13 

3.94±0.12 

1.40 

14656 

h 

2.51±0.08 

0.49±0.01 

3.00±0.08 

3.43±0.08 

1.14 

14657 

k 

2.39±0.13 

0.70±0.02 

3.10±0.13 

4.26±0.12 

1.37 

14657 

h 

3.07±0.10 

0.63±0.01 

3.70±0.10 

4.17±0.10 

1.13 

14658 

k 

1.87±0.13 

0.56±0.02 

2.42±0.13 

3.61±0.12 

1.49 

14658 

h 

2.25±0.08 

0.50±0.01 

2.75±D.08 

3.19±0.08 

1.16 

14659 

k 

2.50±0.13 

0.87±0.02 

3.37±0.13 

4.39±0.12 

1.30 

14659 

h 

3.33±0.08 

0.69±0.01 

4.02±0.08 

4.47±0.08 

l.n 

14672 

k 

1.27±0.15 

0.46±0.02 

1.73±0.15 

2.71±0.14 

1.57 

14672 

h 

1.60±0.09 

0.39±0.01 

1.99±0.09 

2.33±0.08 

1.17 

14673 

k 

2.56±0.14 

0.83±0.02 

3.39±0.14 

4.52±0.13 

1.33 

14673 

h 

2.93±0.08 

0.67±0.01 

3.60±0.08 

4.01±0.08 

1.11 

14674 

k 

2.52±0.13 

0.74±0.02 

3.26±0.13 

4.48±0.12 

1.37 

14674 

h 

2.71±0.08 

0.62±0.01 

3.33±0.08 

3.72±0.07 

1.12 

14708 

k 

3.27±0.13 

1.08±0.02 

4.35±0.13 

5.40±0.11 

1.24 

14708 

h 

3.73±0.09 

0.87±0.01 

4.60±0.09 

4.88±0.08 

1.06 

14709 

k 

3.81±0.13 

0.98±0.02 

4.78±0.14 

5.99±0.13 

1.25 

14709 

h 

4.28±0.12 

0.84±0.02 

5.13±0.12 

5.56±0.11 

1.08 

14710 

k 

2.65±0.10 

0.78±0.02 

3.43±0.10 

4.50±0.09 

1.31 

14710 

h 

2.92±0.05 

0.68±0.01 

3.60±0.06 

3.91±0.05 

1.09 

14713 

k 

2.79±0.09 

0.74±0.01 

3.54±0.10 

4.83±0.09 

1.37 

14713 

h 

2.93±0.06 

0.65±0.01 

3.58±0.06 

3.95±0.06 

1.10 

14715 

k 

3.40±0.09 

0.95±0.01 

4.35±0.09 

5.60±0.09 

1.29 

14715 

h 

3.06±0.06 

0.78±0.01 

3.84±0.07 

4.20±0.06 

1.10 

14716 

k 

3.34±0.07 

0.92±0.01 

4.26±0.07 

5.49±0.07 

1.29 

14716 

h 

3.19±0.05 

0.68±0.01 

3.87±0.06 

4.36±0.05 

1.12 

14717 

k 

3.32±0.10 

0.98±0.02 

4.29±0.10 

5.44±0.09 

1.27 

14717 

h 

3.68±0.06 

0.79±0.01 

4.47±0.06 

4.79±0.05 

1.07 

14718 

k 

2.42±0.13 

0.70±0.02 

3.12±0.13 

4.33±0.12 

1.39 

14718 

h 

2.77±0.06 

0.59±0.01 

3.35±0.06 

3.78±0.06 

1.13 


Table 3: Area- weighted mean cloud Ha absorption EW. 

Values for individual clouds are derived from absorption EW-depth relations. 


Date 

Start UT 

End UT 

Mean Intrinsic Ha EW w 
(channels) (A) 

1988 Dec 16 

12:56 

14:12 

30.76 ± 5.32 

1.82 

± 0.31 

1988 Dec 16 

15:30 

16:13 

45.32 ± 2.11 

2.68 

± 0.16 

1988 Dec 17 

10:35 

11:10 

51.62 ± 3.42 

3.05 

± 0.20 

1988 Dec 17 

13:35 

14:53 

38.56 ± 3.07 

2.28 

± 0.18 

1988 Dec 17 

16:12 

16:44 

52.16 ± 2.29 

3.08 

± 0.14 

1988 Dec 19 

14:48 

16:08 

35.99 ± 1.53 

2.13 

db 0.09 

1988 Dec 19 

17:06 

18:04 

45.50 ± 3.66 

2.69 

± 0.22 

1988 ’Dec 21 

03:37 

04:01 

33.00 ± 0.92 

1.95 

± 0.05 

1988 Dec 21 

04:40 

05:00 

46.39 ± 2.28 

2.74 

± 0.13 

1988 Dec 21 

06:33 

06:46 

42.69 ± 4.74 

2.52 

± 0.28 

1988 Dec 23 

04:53 

05:43 

48.53 ± 2.87 

2.87 

± 0.17 

1988 Dec 23 

05:55 

06:33 

49.14 ± 2.90 

2.90 

db 0.17 

1988 Dec 23 

07:36 

08:01 

50.91 ± 3.01 

3.01 

db 0.18 


Table 4: Atomic parameters and line broadening coefficients 



Ha 

Call K 

Call H 

Mgll k 

Mgll h 

Ae 

- . 

2.27(-6) 

2.27(-6) 

3.34(-5) 

3.34(-5) 

At. wt. 

1.0 

40 

40 

24 

24 

flu 

0.641 

0.66 

0.33 

0.63 

0.31 

A 0 

6562.808 

3933.644 

3968.470 

2795.53 

2802.70 

Tr 

5.698(8) 

1.50(8) 

1.50(8) 

2.70(8) 

2.70(8) 

Gvw 

3.85(-8) 

1.6(-8) 

1.6(-8) 

1.4(-8) 

1.4(-8) 

Gres 

Notes: 

5.877(-7) 


— 

~ 

- 


i- r = r* + Nh(g vw (t/5ooo ) 0 - 3 + c re# ) 


2. Ca II broadening coefficients from Shine & Linsky (1974) 


3. Ha and Mg II broadening coefficients from Lites et ai (1987) 


Table 5: Observed Call H and K index - Ha EW relations 


Date 

UT 

UT 

Echelle 

dI K /dW 

dltf/dW 


Start 

End 

Order 

(A) 


1988 Dec 19 

14:41 

16:19 

143 



0.125 ± 0.027 

1988 Dec 19 

14:41 

16:19 

144 

0.188 ± 0.037 

— 

1988 Dec 19 

16:39 

17:42 

143 



0.071 ± 0.025 

1988 Dec 19 

16:39 

17:42 

144 

0.261 ± cr.071 

— 

1988 Dec 22 

04:12 

04:53 

087 

0.299 ± 0.022 

0.171 ± 0.014 

1988 Dec 22 

04:12 

04:53 

088 

0.256 ± 0.020 

— 

1988 Dec 22 

05:28 

06:03 

087 

0.152 ± 0.021 

0.071 ± 0.013 

1988 Dec 22 

05:28 

06:03 

088 

0.171 ± 0.019 

— 

1988 Dec 22 

07:16 

07:45 

087 

0.300 ± 0.032 

0.207 ± 0.020 

1988 Dec 22 

07:16 

07:45 

088 

0.030 ± 0.029 

— 

1988 Dec 23 

04:47 

05:49 

143 



0.061 ± 0.013 

1988 Dec 23 

04:47 

05:49 

144 

0.083 ± 0.016 


1988 Dec 23 

04:47 

05:49 

145 

0.044 ± 0.023 

— 

1988 Dec 23 

06:01 

06:52 

143 



0.027 ± 0.011 

1988 Dec 23 

06:01 

06:52 

144 

0.119 ± 0.014 

— 

1988 Dec 23 

06:01 

06:52 

145 

0.049 ± 0.020 

— 

1988 Dec 23 

07:17 

07:54 

143 



-0.009 ± 0.021 

1988 Dec 23 

07:17 

07:54 

144 

0.025 ± 0.027 

— 

1988 Dec 23 

07:17 

07:54 

145 

0.050 ± 0.037 

— 



Table 6; Simultaneous Mgll h and k index and Hq EW measures. 

Note that IT is a time-averaged value, calculated over the duration of the corresponding IUE 


Date 



UT 

UT 

W 


Ik 



Ih 





Start 

End 

(A) 







1988 

Dec 

17 

09:51 

10:26 

0.20 

0.00 

± 

0.05 

0.00 

± 

0.03 

1988 

Dec 

17 

11:40 

12:20 

0.57 

0.27 

± 

0.04 

0.29 

± 

0.03 

1988 

Dec 

17 

12:50 

13:30 

0.52 

0.20 

± 

0.04 

0.12 

± 

0.03 

exposure. 19gg 

Dec 

17 

14:04 

14:44 

0.79 

0.38 

± 

0.04 

0.35 

± 

0.02 

1988 

Dec 

17 

16:00 

16:40 

0.49 

0.13 

± 

0.04 

0.04 

± 

0.03 

1988 

Dec 

19 

13:20 

13:55 

0.35 

0.55 

± 

0.04 

0.53 

± 

0.02 

1988 

Dec 

19 

14:29 

15:04 

0.22 

0.13 

± 

0.05 

0.14 

± 

0.03 

1988 

Dec 

19 

15:37 

16:17 

0.43 

0.16 

± 

0.05 

0.21 

± 

0.03 
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Confidence limits for 1988 Dec 17 & 19 Mg II data 





